Self-energy and critical temperature of weakly interacting bosons 
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Using ttie exact renormalization group we calculate the momentum-dependent self-energy at 
zero frequency of weakly interacting bosons at the critical temperature Tc of Bose-Einstein conden- 
sation in dimensions 3 < 7? < 4. We obtain the complete crossover function interpolating between 
the critical regime k <^ kc, where oc k'^~^ , and the short-wavelength regime k 3> fee, where 

E(fc) cx in _D > 3 and cx ln(A;/A;c) in D = 3. Our approach yields the crossover scale kc 

on the same footing with a reasonable estimate for the critical exponent in D = 3. From our T.{k) 
we find for the interaction-induced shift of Tc in three dimensions ATc/Tc « 1.23 an^^'"^, where a is 
the s-wave scattering length and n is the density. 
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Introduction. Triggered by Nobel Prize winning exper- 
iments, the past few years have seen an upsurge of ac- 
tivity in the field of Bose-Einstein condensation Q . The 
interacting Bose gas falls into the universality class of the 
classical 0{2) model. While for dimensions D slightly be- 
low the upper critical dimension D = A the small param- 
eter e = 4 — D can be used to study the critical behavior 
via controlled perturbative renormalization group (RG) 
calculations no such small parameter exists in the 
physically relevant case D = 3. The problem of Bose- 
Einstein condensation in a three-dimensional Bose gas 
therefore requires non-perturbative methods. 

One of the fundamental properties of the system is the 
energy dispersion of the elementary excitations. In 
the vicinity of the critical temperature Tc, where the cor- 
relation length becomes large, interactions qualitatively 
modify the energy dispersion at long wavelengths. Pre- 
cisely at Tc it is well known that ek oc fc^~'' for sufficiently 
small momenta k. In D = 3 high order re-summations 
of the e-expansion 3] and recent Monte-Carlo simula- 
tions 4] yield the critical exponent rj w 0.038. On the 
other hand, when k exceeds some interaction-dependent 
crossover scale fee, self-energy corrections to the non- 
interacting dispersion p^k^ become negligible. Here po — 
h? /2m, with some effective mass m. Let us describe this 
crossover in terms of a dimensionless self-energy function 



a{x) = (pofcc)"MS(fcca:) - S(0)] 



(1) 



where S(fc) is the exact self-energy of the interacting Bose 
gas at zero frequency. By construction, the crossover 
occurs at x — 0{1). Note that the field-theoretical RG 
effectively sets fee oo, so that neither kc nor cr(x) for 
a; 3> 1 can be obtained within that framework. 

Naively one might expect that in the short-wavelength 
regime x ^ 1 the scaling function can be calculated per- 
turbatively. Unfortunately, in Z? = 3 one encounters in- 
frared divergences even in this regime. Baym et al. 5] 
recently employed a simple self-consistent re-summation 



which cures the divergences in _D = 3. However, their 
result 77 = 1/2 is far from the above mentioned accepted 
value 13, 01 ■ Clearly, it would be very useful to have 
a quantitatively accurate formula for the momentum- 
dependence of the self-energy at Tc, which interpolates 
between the critical and the short-wavelength regime. As 
emphasized in Ref. Q, the shift in the critical tempera- 
ture ATc due to interactions is closely related to the be- 
havior of the scaling function cr(x); at weak coupling and 
to leading order in the scattering length, ATc in _D = 3 
is proportional to the crossover scale kc, with a constant 
of proportionality that depends on the specific form of 
a{x) for < X < oo. 

Exact RG flow equations for interacting bosons. Previ- 
ously Bijlsma and Stoof Q used RG methods to investi- 
gate the infrared behavior of the weakly interacting Bose 
gas above and below Tc. However, these authors consid- 
ered only the constant parts of the vertices. In particu- 
lar, no attempt was made to calculate the momentum- 
dependence of the self-energy. Here we study this prob- 
lem using the exact RG in the form proposed by Wet- 
terich Q, and by Morris 8]. In this approach one con- 
siders the RG flow of the generating functional of the 
one-particle irreducible vertices F^^"-* as an infrared cut- 
off A is reduced. Assuming that the bosons are in the 
normal state 9], the exact flow equation for the irre- 
ducible two-point vertex is [lol lll| 



9Arf (if) 



K' 



where for sharp infrared cutoff A 

Si\k\-A) 



Ga{K) = 



zc^„-ek + M-S(0,zO)-rf (K) 



(2) 



(3) 



We use the notation K = (k, zwn), — {(iV)^^ ^ , 
where (3 is the inverse temperature, V is the volume, 
and ujn — 27toT are bosonic Matsubara frequencies. The 
irreducible four-point vertex satisfies 
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GK{K)Tf\K[,K'^, K; K, K2,Ki) 
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(4)/ 



K'=Ki+K2-K 



[G^{K)GdK') + GA{K)GA{K')\Tf(K[,K'; K, Ki)T^^\k^, K- K' , K2) 
[Ga{K)Ga{K') + Ga{K)Ga{K')]t'^^\k^,K'; K, Ki)T^^\k[, K; K', K2) 

I 



K'=Ki-K{+K 



K'=Kx-K'+K 



(4) 



Here the exact propagator Ga{K) is obtained from Eq. 
© by replacing (5(|k| - A) ^ e(|k| - A). The flow equa- 
tion for the irreducible six-point vertex rj^' on the right- 
hand side of Eq. Q can be found in Ref. Be- 
cause quantum mechanics is irrelevant for classical crit- 
ical phenomena, for our purpose it is sufficient to retain 
only the zero frequency part of all vertices and work with 
an effective classical theory, with ultraviolet cutoff given 
by Ag — 27r/Ath, where Ath = h / 2innT is the ther- 
mal de Broglie wavelength. In principle, the above flow 
equations can also be used to calculate the parameters 
of the effective classical model, which are renormalized 
by the degrees of freedom involving non zero Matsubara 
frequencies. Here we take these finite renormalizations 
implicitly into account via the initial conditions at scale 
Aq. In particular, in = 3 we write the initial value of 
the four-point vertex as (0,0; 0,0) = I&tt p^a, where 
a is the s-wave scattering length. 

To study the critical regime, it is convenient to intro- 
duce dimensionless momenta q = k/A and the rescaled 
classical propagator Gi{q) = — (poA^/Zi)G'A(qA, zO) = 
Q{q—1)/ Ri (q) , which is a function of the logarithmic flow 
parameter I = -ln(A/Ao). Here Ri{q) = Ziq^ +ff\q), 



where tf\q) = Zi{p„K^) 

two-point vertex, and Zi = 1 
wave-function renormalization factor. For n > 2 we de- 
fine the rescaled higher-order vertices 



^-lp(2) 

^ A 



(Aq, iO) is the rescaled 
21 



(.<l)/dq\2=o is the 



fP")(q,) = (if,3r)"-iA^("-i)-2" (Z,/p„)^ 



X F^^"'(k, -> Aqj,u;„ 



0) 



(5) 



For later convenience we have inserted the numerical fac- 
tor Ko = ri£)/(27r)-^, where is the surface area of the 
D-dimensional unit sphere. The rescaled vertices satisfy 
exact flow equations of the form 



2n — D{n — 1) — nrji 



q', • + q. • Vq.) f [^")(qO + rp"^(qO , (6) 



1=1 



p(2n). 



where rji = dr'i'^\q)/dq'^\q=o is the flowing anomalous 
dimension. The function F; (g) on the right-hand side 



of the flow function of the two-point vertex is 



rp'(.) 



Gi{q')T[^\q,q';q',q) 



(7) 



with = / d^q/rtD and Gi{q) = 6{q ~ l)/Ri{q). The 

flow function f P'' (q'l, q2; q2, Qi) can be obtained from 
the right-hand side of Eq. @ by replacing Ga{K) 

-Gi{q). Ga{K) ^ -Gi{q). rP"^ - f P"\ 4, and 

multiplying the resulting expression by an overall minus 
sign. From the above deflnitions we can show that 
at Tc the function Q can be written as 



a{x) 



d/e"2(i-ic)+/o'<ir,Kf(2- 



(8) 



where F 



tf\0) and Ic 



will be deter- 



mined shortly. Hence, the flow function Q is the key 
quantity for obtaining a{x) via the RG. 

Classification of couplings and crossover scale. In 3 < 
D < A there are two relevant couplings: the momentum- 
independent parts ri = ff\Q) and ui = f^^ (0,0; 0,0) of 
the two-point and four-point vertices. If we simply ignore 
all other couplings, the flow of r; and ui is determined by 



dm = {2-rii)ri+ui/{l + ri) , (9) 
dm ^ {€~2Tji)ui~'^uy{l + rif . (10) 

Note that the initial value Tq has to be fine tuned so 
that r^ = limj^oo fi is finite, corresponding to the criti- 
cal point. As long as |r/| ^ 1 and 77/ ^ e, the qualitative 
behavior of the solution of Eq. H10() can be obtained 
by simply replacing r; ^ and 77; — > on the right- 
hand side. Then the solution of Eq. 1)10(1 is a Fermi 
function, = [e<^('<:~') + 1]^^7 where = 2e/5 and 

Ic = hi{uif / uq — 1) « e^^ ln(u*/uo). Throughout 
this work we assume that uq ^ u*, corresponding to 
weak bare interactions. Thus, within a narrow interval 
of width 1/e centered dX I — Ic the coupling ui suddenly 
approaches its fixed point value u■^,. Below we shall show 
that kc = Aoe~''= ~ Ao(uo/w*)^^'^ is the scale where the 
dispersion crosses over from k'^ to the k'^~^. 



3 



We now discuss the marginal couplings. For 3 < D < 4 
the only marginal parameter is the wave-function renor- 
malization Z;; it satisfies by construction diZi = —rjiZi. 
In D = 3, there are three additional marginal couplings. 
The most obvious one |^ is the momentum-independent 
part of the six-point vertex, vi = f [^^(0, 0, 0; 0, 0, 0). But 
there are two more marginal couplings ai and bi related to 
the four-point vertex, which has the expansion for small 
momenta 



f[''^(q'i>q2;q2,qi) 



ui + ai (|qi 

^/|qi + q2| - 



-qil + |qi 

0(q?) . 



q2l 



(11) 



One easily verifies that in D = 3 the scaling dimensions 
of ai and bi vanish. Note that even if initially f |1q — uq 
is momentum-independent (corresponding to the initial 
conditions oq = &o — 0), finite values of ai and 6; are 
generated as we iterate the RG. 

Truncation and anomalous dimension. To calculate 
rji, we need the momentum-dependent part of the four- 
point vertex. Note that in I? > 3 the momentum- 
dependence of f 1^ Vq'i iC(n ; q2 , qi ) is irrelevant. As ex- 
plained in Refs. [ij, [ij], for sufficiently large I irrel- 
evant couplings become local functions of the relevant 
and marginal ones. Assuming u/ <C 1, the momentum- 
dependence of f I**-* {q'l , q2 ; q2 J qi ) can be obtained by set- 



ting r 



(6) 



andr 



(4) 



ui in the expression for the flow 



function f ['^■'(q'^, q'j; q2, qi). In Ref. we have used 
a similar truncation to calculate the spectral function of 
the Tomonaga-Luttinger model using the exact RG. In 
the present problem this truncation amounts to 



rr^(q'i,q2;q2,qi) 



^x/(|qi 



+Xi(|qi - qil) + Xi(|qi 



+ q2|) 
q2l)" , 



(12) 



where the generalized susceptibility is given by xi{q) = 
2 / ,Gi{\c['\)Gi{\c[' + q|). The approximation ((EJ can 
be formally justified as long as ui remains small, which is 
certainly the case for small e. Unfortunately, in I? = 3 the 
renormalized m, is not guaranteed to be small. Moreover, 
the couplings vi, ai and bi which become marginal in 
= 3 are not consistently taken into account in Eq. 
()12|l . so that a priori we cannot estimate the accuracy of 
our truncation in D = 3. Nevertheless, we shall use Eq. 
(|12|l as a first approximation even in Z) = 3; an improved 
truncation including the marginal couplings vi , ai and bi 
will be presented elsewhere (l^ . 

Given the truncation H12|l. the exact flow equation © 

for f l^"* (q'l, q2; q2, qi) can be solved without further ap- 
proximation. Substituting the result into Eq. Q we 
obtain for the fimction rp^'^'(g) required in Eq. © 



(9) 



x(x;_,(e-*|q' + q|)-x,_i(e-*))q' , (13) 



where (. . .)q' denotes angular average over the unit vector 
q'. Eq. (|13() implies the following integral equation for 
the flowing anomalous dimension rji = dt'f^ [q) / dq^ ' 



1^1 = / dtK(l,t)u'i_^e-^^i-^'''''- 
Jo 



(14) 



with the kernel given by 
3 



Kil,t) 



(15) 



where x'i{x) = dxi(x)jdx and j(!l{x~) = d"^ xi{x) / dx"^ . Eq. 
(|14|l together with the two flow equations (|9I10() form a 
system of three coupled integro-diffcrcntial equations for 
the three unknown functions n, ui and 77/. For small e 
we recover the known result 77 — lim/_»oo '?/ = e^/50. In 
Z) = 3 we find numerically at the fixed point r.^, — —0.143, 
— 0.232, and = 0.104, which should be compared 
with the accepted value 77 = 0.038 |^ Given the 

simplicity of our truncation, it is quite satisfactory that 
our estimate for 77 has the correct order of magnitude. 
We can substantially improve our estimate for 7? by in- 
cluding the above mentioned marginal couplings [iq . In 
comparison, the method of Ref. Q yields r] — 0.5. 

Scaling function. To obtain the scaling function a{x) 
we substitute our approximation p3|l into Eq.®. The 
integrations can be simplified by noting that the flow of 
77r on the right-hand side can be approximated by a Fermi 
function, 77^. = 7;[e'''^~'^'' -I- 1]~^. We have verified numer- 
ically that an analogous substitution in Eq. (fT^ yields a 
good approximation to the true solution rji . Introducing 
appropriate integration variables we obtain 



a{x) 



3ui 



dy 



y ^+^'F{x,y;Tj,Q 

[x" + ?7e]2-277/£ 



(16) 



with the dimensionless function 



F{x,y;'q,lc) 



dz- 



"'[x^ + {y/zYf'[l + ri^+,,,^] 

X {Xh+\n I {z) - X/c+ln I (kq' + 2/q|)>q' • (17) 

A numerical evaluation of Eq. H16|) is shown in Fig. ^ 
The asymptotic behavior for small and large x can be 
extracted analytically. In the critical regime x <^ 1 
we find in D dimensions (j{x) ~ A^x"^^"^ , with Ajj = 
^ulj^ dyy~^+'^'^F{Q,y;ri,oo). In Z) = 3 we obtain 
numerically A3 — 1.20. On the other hand, in the 
short wavelength regime we find cr{x) - Bdx^'-"-^^ for 
[2{D - 3)]"^ < Inx, and a{x) ^ iJglna; for 1 < Inx < 
[2{D -3)]-^ (which includes the limit D ^ 3). Here 
Bd oc ul/{D - 3) and B'^ = 3Tr'^ul/2A. Keeping in mind 
that kc cx (uo/utY^'^, it is easy to see that for k ^ kc the 
physical self-energy E(fc) — T,{0) is proportional to Uq. 
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f 




FIG. 1: Double logarithmic plot of the scaling function a{x) 
in D = 3, see Eqs. ^ and H16|l . Note that in our rescaled 
units we obtain a universal curve for uo — > 0. The dashed 
lines indicate the asymptotic behavior for small and large x. 

Critical temperature. Given a'(x), we may calculate the 
interaction-induced shift ATc of the critical temperature 
of the weakly interacting Bose gas. Generalizing an ex- 
pression given by Baym et al. Q for arbitrary D, we 
obtain for the contribution from classical fluctuations 

T, 7rDC{D/2) Jo + a{x) " ^ ' 

Keeping in mind that a{x) ~ Bdx-^^^~^^ for Z? > 3 and 
large a:, it is easy to see that we may remove the ultra- 
violet cutoflF (Aq oo) in Eq.l^HJ) provided |^ < 2, 
i.e. D < 10/3. Only in this case AT, is dominated by 
classical fluctuations. For D > 10/3 the value of the 
integral in Eq. H18|) depends on the ultraviolet cutoff, 
such that it is proportional to Uq, with logarithmic cor- 
rections (cx UqIumo) in D = 10/3. In three dimensions, 
where uo = 167r-i[C(3/2)]-i/3a„i/3^ we find [T^ to lead- 
ing order ATc ~ cian^^^Tc with ci = 1.23, in agreement 
with the five-loop result ci = 1.14±0.11 by Kleinert [l^ 
and the seven-loop result ci = 1.27 ± 0.11 by Kastening 
[l^. For recent reviews on the problem of calculating 
AT, see Ref. H^. 

Summary and outlook. In this work we have used the 
exact RG to calculate the momentum-dependent self- 
energy of weakly interacting bosons at the critical point. 
Our central result is an explicit expression for the scaling 
function a{k/kc) (see Eqs. ^ and which interpo- 

lates between the critical regime k ^ kc and the short- 
wavelength regime k ^ kc. Note that the momentum- 
dependent part of Eq. H12|l contains an infinite set of 
irrelevant (and hence non-rcnormalizable) couplings, so 
that neither the crossover scale kc nor S(fc) for k > kc 
are accessible within field-theoretical RG. The crossover 
function shown in Fig.n]should be useful for a comparison 
with experiments probing the single-particle excitations 
in a wide range of length scales, for example via mea- 



surements of the velocity distribution. We have shown 
here how the exact RG can be employed to calculate full 
scaling functions and not just their asymptotic behavior. 
Our method is quite general and should prove useful in 
the study of other critical systems, for example Luttinger 
liquids [l3| or systems described by the dilute Bose gas 
quantum critical point |2l|. 

This work was partially supported by the DFG via 
Forschergruppe FOR 412. 
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